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ABSTRACT 

We study the relation between accretion, black hole mass and jet power in AON, by using a large group of 
blazars detected by the Eermi Large Area Telescope and radio galaxies. Our main results are as follows, (i) The 
jet power of ESRQs and ERII-HEG depends on the black hole mass, which suggests that the ESRQs and ERII- 
HEG are in Radiation-Pressure Dominated regime. The jet power of BL Lacs and ERI-LEG depends on the 
accretion, which suggests that the BL Lacs and ERI-LEG are in the Gas-Pressure Dominated regime, (ii) We 
hnd that most of ESRQs and BL Lacs have Pjet > Lg which suggests that the Blandford-Znajek mechanism 
is insufficient to explain the jet power of these objects, (iii) The ESRQs are roughly separated from BL Lacs 
by the Ledlow-Owen’s dividing line in the log Pjet — log M plane, which supports the unified scheme of AGN. 
(iv) The ESRQs and BL Lacs have a clear division at Lboi/LEdd ~ 0.01, and the distribution of Eddington 
ratios of BL Lacs and ESRQs exhibits a bimodal nature, which imply that the accretion mode of ESRQs may 
be different from that of BL Lacs, (v) We find a significant correlation between broad line luminosity and jet 
power, which supports a direct tight connection between jet power and accretion. 

Subject headings: galaxies: jets-galaxies: active-BL Lacertae objects: general-quasars: general- 
accretion,accretion disk 


1. INTRODUCTION 

Blazars are the most extreme Active Galactic Nuclei (AGN) 
with a relativistic jet directed towards our line of sight (Urry & 
Padovani 1995). The emission from the jets is highly boosted 
and dominates the AGN emission at all wavelengths due to 
a relativistic beaming effect (Sbarrato et al. 2014). Blazars 
are generally divided into two subcategories including BL Lac 
objects (BL Lacs) and Plat spectrum radio quasars (ESRQs). 
BL Lacs only have very weak or no emission lines, whereas 
ESRQs show strong emission lines. The classical division be¬ 
tween ESRQs and BL Lacs is mainly based on the rest frame 
equivalent width (EW) of their broad emission lines. Partic¬ 
ularly, if the EW is larger than SA, the objects are classified 
as ESRQs, otherwise as BL Lacs (Urry & Padovani 1995). 
The broad emission lines of ESRQs are produced by distant 
gas clouds in broad-line regions (BLRs) which are photoion- 
ized by the optical/UV continua radiated from the accretion 
disks surrounding massive black holes (Sbarrato et al. 2012). 
The difference of the broad line emission between ESRQs 
and BL Lacs may be attributed to their different central en¬ 
gines (Cavaliere & D’Elia 2002; Cao 2002, 2003). Landt et 
al. (2004) introduced an analogous classification criterion and 
found that it is possible to discriminate between objects with 
intrinsically weak or strong narrow emission lines by study¬ 
ing the [On] and [Om] EW plane. Abdo et al. (2010a) 
have classified blazars based on the synchrotron-peak fre¬ 
quency of the broadband SED. Ackermann et al. (2011) used 
the estimated value of Vp^^j. to classify the sources as ei¬ 
ther a low-synchrotron-peaked blazar (LSP, for sources with 
Vpeak < 10 '^"'Hz), an intermediate-synchrotron-peaked blazar 
(ISP, for lO^^Hz < v^aak < lO^^Hz), or a HSP blazar (if 
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V^eak > 10'^ Hz). 

Radio galaxies are often thought to be the parent popula¬ 
tion of blazars. The blazars are aligned to our line of sight, 
radio galaxies have their jets oriented at larger viewing an¬ 
gle (Urry & Padovani 1995; Sbarrato et al. 2014). Radio 
galaxies are commonly divided into two subclasses includ¬ 
ing PRI and PRII based on their radio morphology (Panaroff 
& Riley 1974). PRI shows bright jets close to the nucleus, 
while PRII shows prominent hot spots are far from it. The 
PRI and PRII radio galaxies can be clearly divided in the 
host galaxy optical luminosity-radio luminosity (Mr — Lj-ad) 
plane by a dividing line showing that radio power is propor¬ 
tional to the optical luminosity of the host galaxy (Ledlow 
& Owen 1996). Wu and Cao (2008) found that the dividing 
line of the Ledlow-Owen relation for PRI/PRII can be roughly 
reproduced by using the maximum jet power available from 
ADAPs around kerr black hole as a function of black hole 
mass with a hybrid jet formation model (i.e., BPh-BZ mech¬ 
anism). However, what causes the PRI/PRII division is still 
unclear. Gopal-Krishna and Wiita (2000) suggested that the 
differences between PRI and PRII may be caused by the 
jets interact with the ambient medium with different physi¬ 
cal properties. However, some authors thought that the differ¬ 
ences between PRI and PRII may be the intrinsic difference 
of their central engines, such as different accretion models 
and the formation process of jets (e.g., Baum et al. 1995; 
Bicknell 1995; Reynolds et al. 1996; Meier 1999; Ghisellini 
& Celotti 2001; Marchesini et al. 2004; Hardcastle et al. 
2007). Ghisellini and Celotti (2001) found that the division 
between PRI and PRII actually reflected a systematic differ¬ 
ence in accretion rate. The PRIs have generally lower accre¬ 
tion rate (rh < 0.01), while for PRII it was typically larger. 
Some authors also found that the ESRQs and BL Lacs can 
be divided in rh ~ 0.01 (see e.g. Narayan, Carcia & McClin- 
tock 1997; Ghisellini, Maraschi & Tavecchio 2009). They 
suggested that the PRI and PRII could correspond to differ¬ 
ent accretion models. Laing et al. (1994) suggested that 
PRIIs can be divided into high-excitation (HEG) and low- 


2 


excitation galaxies (LEGs), which are also applied to some 
FRI. Buttiglione et al.(2009,2010) found that all HEG are 
ERIE while LEGs can be both FRI and ERIE Best & Heck¬ 
man (2012) suggested that HEG radio galaxies are fuelled 
at high rates through radiatively-efficient standard accretion 
disks, while LEG radio galaxies are fulled via radiatively in¬ 
efficient flows at low accretion rates. Heckman & Best (2014) 
suggested that most FRIs are LEG sources, and most FRIIs 
are HEG sources. Gendre et ak (2013) suggested that the 
LEG and HEG sources could be divided by the Ledlow and 
Owen’s line. These results may suggested that the FSRQs vs 
FRII and BE Lac vs FRI unification can be easily reinterpreted 
physically as HEG/LEG because of the large population over¬ 
laps between HEGs and FRII and between LEGs and ERE 

In the frame of unification schemes, many authors have 
studied the unified scheme of FSRQs and FRII radio galax¬ 
ies. Padovani and Urry (1992) found that the radio luminosity 
functions (RLFs) of FRII radio galaxies are consistent with 
the observed radio luminosity functions of FSRQs, which fa¬ 
vor the unification of FRII radio galaxies and FSRQs. Anal¬ 
ogously, the unified scheme of BE Lac objects and FRI radio 
galaxies has been extensively studied by many authors who 
used different approaches, such as comparison of spectral en¬ 
ergy distributions (SEDs) in different wavebands and radio 
morphology (Owen et al. 1996; Capetti et al. 2000; Bai & 
Lee 2001). Padovani and Urry (1991) derived that the RLFs 
of BE Lacs are consistent with the observed RLFs of FRI ra¬ 
dio galaxies, which also favor the unification of BE Lacs and 
FRI radio galaxies. Since the launch of the Fermi satellite, 
we have entered in a new era of blazars research (Abdo et 
al. 2009, 2010). Up to now, the Fermi Large Area Telescope 
(EAT) has detected hundreds of blazars because it has about 
20 folds better sensitivity than its predecessor EGRET in the 
0.1-100 GeV energy rang. Sbarrato et al. (2014) found a good 
correlation between the broad line luminosity and radio lumi¬ 
nosity, which suggested a direct tight connection between jet 
power and accretion rate. They got that the observational dif¬ 
ferences between blazar subclasses reflect differences in the 
accretion regime. 

Blandford and Rees (1978) suggested that the radiation ob¬ 
served from blazars is dominated by the emission from rela- 
tivistiv jets, which transport energy and momentum to large 
scales. However, the formation of jets is not still unclear 
in astrophysics (e.g. Meier et al. 2001). Many models 
have been proposed to explain the origin of jets. There are 
two main theoretical models of the formation of jet. One 
is that the rotational energy of the black hole is expected to 
be transferred to the jets by the magnetic fields threading the 
holes when assuming the black hole is spinning rapidly (i.e., 
the Blandford-Znajek (BZ) mechanism; Blandford & Znajek 
1977). Another one is that the jet can also be accelerated by 
the large-scale fields threading the rotating accretion disk (i.e., 
the Blandford-Payne (BP) mechanism; Blandford & Payne 
1982). In the BZ theory. Gosh & Marek (1997) found that 
it is possible to find two regimes for standard disc, which 
may be dominated by the radiation pressure (RPD) or the gas 
pressure (GPD). Foschini (2011) suggested that the FSRQs 
and narrow-line Seyfert 1 are in RPD regime, while BE Lacs 
are basically in the GPD regime by using the data of 30 FS¬ 
RQs and 9 BE Lacs from Ghisellini et al. (2010) and 4 7 -ray 
narrow-line Seyfert 1 from Abdo et al. (2009a). 

In this paper, we use a sample of blazars detected by Fermi 
Large Telescope (EAT) and radio galaxies with measured jet 
kinetic power, black hole mass and Eddington ratio to explore 


the relationship between them. The samples are described in 
Section 2; the results and discussions are presented in Section 
3; conclusions are in Section 4. The cosmological parameters 
Hq = 70kms”^Mpc“^, Hm = 0.3, and Ha = 0.7 have been 
adopted in this work. 

2. THE SAMPLE 

In order to study the relationship between FSRQs, BE Lacs, 
FRI and FRII radio galaxies, we tried to select the large group 
of clean Fermi blazars and radio galaxies with reliable red- 
shift, black hole mass and jet kinetic power. 

2 . 1 . The Fermi blazars 

We tried to get the large group of Fermi blazars with reliable 
broad line luminosity, black hole mass and jet kinetic power. 
We consider the sample of Xiong & Zhang (2014). Xiong 
and Zhang (2014) collected a larger group of Fermi blazars 
with broad line luminosity, black hole mass and jet power. 
However, there are some blazars that we don’t know the types 
in their samples. Hence, we only collected known types of 
blazars from their samples. Firstly, Xiong and Zhang (2014) 
considered the following samples to get broad line data: Cao 
& Jiang (1999), Wang et al. (2004), Liu et al. (2006), Sbar¬ 
rato et al. (2012), Chai et al. (2012), Shen et al. (2011), 
Shaw et al. (2012). Secondly, they considered the following 
samples to get black hole mass: Woo & Urry (2002), Xie et 
al. (2004), Liu et al. (2006), Zhou & Cao (2009), Zhang et 
al. (2012), Sbarrato et al. (2012), Chai et al. (2012), Leon- 
Tavares et al. (2011a), Shen et al. (2011), Shaw et al. (2012). 
Nemmen et al. (2012) estimated the jet power for a large sam¬ 
ple of Fermi blazars. Xiong and Zhang (2014) collected the 
jet power from the Nemmen et al. (2012) for their sample. 
When the jet powers of blazars were not directly gotten from 
the Nemmen et al. (2012), Xiong and Zhang (2014) used the 
relationship between intrinsic 7 -ray luminosity and jet power 
derived by Nemmen et al. (2012) to calculate the jet power. 
We use the broad line luminosity to get bolometric luminos¬ 
ity (Lboi ~ IOLblr) and use the ratio of bolometric luminos¬ 
ity and Eddington luminosity to calculate the Eddington ra¬ 
tio (rh = Lboi/LEdd, LEdd = 1.3 X 103®(M/MQ)ergs ^ for 
Fermi blazars. At last, we get a sample containing 159 clean 
Fermi blazars (112 FSRQs and 47 BE Lacs). 

2.1.1. The broad line luminosity 

Celotti et al. (1997) calculated the broad-line luminosity by 
scaling several strong emission lines to the quasar template 
spectrum of Francis et al. (1991), using Lya as a reference. 
Francis et al. (1991) set Lya flux contribution to 100, and the 
relative weights of the Ha, H/3, Mgll and CIV lines to 77, 
22, 34, and 63, respectively. The total broad line flux is fixed 
at 555.76. The broad line luminosity value of each source 
has been derived by using these properties. When more than 
one line was presented, the approach that calculated simple 
average of broad-line luminosity estimated from each line was 
adopted. The rest of authors in our sample also used method 
proposed by Celotti et al. (1997) to calculate the broad line 
luminosity. 

2.1.2. The black hole mass 

Usually, there are three methods to calculate black hole 
mass, namely, the traditional virial black hole mass, the stellar 
velocity dispersion or the bulge luminosity, and the variation 
timescale. The black hole mass is estimated by traditional 
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virial method for most of FSRQs in the sample (Woo & Urry 
2002; wang et al. 2004; Liu et al. 2006; Sbarrato et al. 2012; 
Chai et al. 2012; Shaw et al. 2012; Shen et al. 2011; Wu et al. 
2011). For some sources, especially BL Lacs, the black hole 
mass can be estimated by the stellar velocity dispersion or 
the bulge luminosity (Woo & Urry 2002; Zhou & Cao 2009; 
Zhang et al. 2012; Sbarrato et al. 2012; Chai et al. 2012; 
Leon-Tavares et al. 2011a). For few sources, the black hole 
masses are estimated by the variation timescales (Xie et al. 
1991, 2004, 2007). When more than one black mass is gotten 
for same objects in the sample, the average black hole masses 
are adopted. 

We also notice the possible biases due to different meth¬ 
ods were used for different samples. The black hole mass was 
calculated by different methods in our sample. Tremaine et al. 
(2002) suggested that the uncertainty in the M — cr relation is 
small, <0.21 dex. The uncertainty on the zero point of the 
line width-luminosity-mass relation is approximately 0.5 dex 
(Gebhardt et al. 2000; Ferrarese et al. 2001). McLure & Dun¬ 
lop ( 2001 ) suggested that the uncertainty of the Mbh — Mr 
relation is 0.6 dex. We assume that the average uncertainty of 
black hole mass is 0.44 dex. The estimated black hole masses 
from stellar dispersion velocity are log(M/M0)=8.29 (Mkn 
421), 9.21 (Mkn 501), 9.03 (2005-489), 8.23 (BL Lac) (Woo 
& Urry 2002). The traditional virial black hole masses are 
log(M/MQ)=8.5 (Mkn 421), 9 (Mkn 501), 8.5 (2005-489), 
8.7 (BL Lac) (Sbarrato et al. 2012). The estimated black 
hole mass from the variation timescales are log(M/MQ)=7.6 
(Mkn 421), 8.3 (Mkn 501), 8.1 (2005-489), 7.8 (BLLac)(Xie 
et al.2004). We find that the largest difference between them 
is less than 1 order of magnitude. Furthermore, we also find 
that the black hole masses were calculated by the variation 
timescales are a little lower than the black hole mass calcu¬ 
lated by another two methods. However, these sources are 
very few, and have no significant impact on our results. 

2.1.3. The jet kinetic power 

The jet kinetic power is estimated by using the correlation 
between the extended radio emission and the jet power. Ac¬ 
cording to search for X-ray cavities in different systems in¬ 
cluding giant elliptical galaxies and cD galaxies, Cavagnolo 
et al. ( 2010 ) obtained this tight correlation between the “cav¬ 
ity” power and the radio luminosity 

P- “ (1) 

which is continuous over ^ 6-8 decades in Pjet and Pradio 
with a scatter of «0.7 dex and Pjet = Pcav Meyer et al. 
( 2011 ) used this formula to calculate the jet power for ra¬ 
dio galaxies. Following Meyer et al. (2011), Nemmen et al. 
( 2012 ) also adopted this formula to calculate the jet power for 
a large sample of Fermi blazars. They got the relationship 
between the intrinsic 7 -ray luminosity and the jet power 

log Pjet = (0.98 ± 0.02) logL“*-f (1.6 ±0.9) ( 2 ) 

with a scatter of ^0.64 dex. Nemmen et al. 

( 2012 ) also obtained the relationship between ob¬ 
servation 7 -ray luminosity and beaming factor, 
fb ~ 5 X 10“‘^(Lobs/10'^®ergs“^)“°-^®^°'^®, where fb 

is a beaming factor. Lobs is observation 7 -ray luminosity. The 
intrinsic 7 -ray luminosity L can be calculated by L = fbLobs- 
When the jet powers of blazars were not directly gotten from 
the Nemmen et al. (2012), Xiong & Zhang (2014) used the 
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Figure 1. The jet power calculated by equation (1) as a function of jet power 
calculated by equation (2) for FSRQs (black), LBL (red), HBL (green), IBL 
(blue). The black line is log pgamma 

equation (2) to calculate the jet power. We should note the 
low significant of the fb vs Lobs relation derived by Nemmen 
et al. ( 2012 ). 

We also note that the jet power was calculated by two dif¬ 
ferent methods, namely, the intrinsic radio luminosity (using 
a relationship with 0.7 dex scatter) and the intrinsic 7 — ray 
luminosity (0.6 dex). We compare the jet power that calcu¬ 
lated by equation ( 1 ) with that calculated by the equation ( 2 ) 
in Figs.l. Figure 2 shows that the jet power that calculated 
by equation ( 1 ) and that calculated by the equation ( 2 ) as a 
function of black hole mass. From figure 2, we can see that 
the distributions of jet power and black hole mass have no sig¬ 
nificant difference for our sample although the jet power was 
calculated by equation (2). From these figures, we can see that 
there is no bias in our results obtained within 3cr. Moreover, 
the bulk of our samples (about 60%) have jet powers that cal¬ 
culated by equation ( 1 ) and the minority of the sample (about 
40%) have jet power that calculated by equation (2). 

2.2. The radio galaxies 

We tried to select the FRll-HEG and FRl-LEG radio galax¬ 
ies with reliable redshift, black hole mass and jet power. 
Eirstly, we consider the sample of Meyer et al. (2011). 
Meyer et al. (2011) estimated the jet power by the method 
of Cavagnolo et al. (2010) for a sample of radio galaxies. 
Buttiglione et al. (2010) have classihed the sources as HEG, 
LEGs and broad line objects (BLOs) according to optical fea¬ 
tures. We cross-correlated Meyer’s sample with Buttiglione’s 
ample. We get the 14 ERl/LEG and 6 ERll/HEG radio galax¬ 
ies. Secondly, we use the Mbh — Mr relation to calculate 
the black hole mass by using an average color correction of 
R—H=2.5 (Dunlop et al. 2003). At last, the narrow-line re¬ 
gions (NLRs) are believed to be photoionized by the radi¬ 
ation from the accretion disk and the narrow-line emission 
can be used to estimate the bolometric luminosity for ra¬ 
dio galaxies. We calculate bolometric luminosity by using 
the relation proposed by Heckman & Best (2014), namely, 
Lboi = 600 L[ori]. We collect the narrow lines [Om] from 
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Figure 2. The jet power calculated by equation (1) and jet power calculated 
by equation (2) as a function of black hole mass for FSRQs (black), LBL 
(red), HBL (green), IBL (blue). 

Buttiglione’s ample and calculate Eddington ratio. 


3. RESULTS AND DISCUSSIONS 


In the BZ theory, the magnetic field of the disc is pushed 
toward the event horizon by the Maxwell pressure (Foschini 
2011). For a standard disc, its magnetic field depends on the 
accretion rate and it is possible to find two regimes (Ghosh 
& Abramowicz 1997). One refers to strong accretion disc, 
which is dominated by the radiation pressure (RPD). Ghosh 
and Abramowicz (1997) pointed out that the jet luminosity of 
the BZ mechanism in the radiation pressure dominated can be 
estimated by the following formula 

Lbz.rpd = 2 X 10'^'^( ^Q8]^^ )j^ergs~^ (3) 

where M is the black hole mass, and j is the dimensionless 
angular momentum. The other is the standard disc with low 
accretion dominated by the gas pressure (Ghosh & Abramow¬ 
icz 1997), which can be expressed as 


LbZ.GPD = 8 X 1044(_^)ll/10(^)4/5j2^^g3-l 


(4) 


where rii is the accretion rate. 


3.1. Jet power versus black hole mass and Eddington ratio 

Figure 3 shows jet power as a function of black hole mass 
(left panel) and jet power as a function of Eddington ratio 


(right panel). From Fig.3, we can see two different regimes. 
One is that jet power depends on the black hole mass and 
another is that jet power depends on the accretion. We find 
that the FSRQs and FRII-HEG are in the mass-dependent 
regime, while BE Facs and FRI-FEG are in the accretion- 
dependent regime. In Fig.3, it is easy to recognize that BE 
Facs and FRI-FEG are basically in the GPD regime, while 
FSRQs and FRII-HEG are in the RPD regime. Observational 
evidence that quasars can accelerate high-velocity winds is 
plentiful, Cattaneo et al.(2009) suggested that this winds is 
‘momentum-drive’ by radiation pressure. According to the 
SED modeling, Foschini (2011) also found that the FSRQs 
and 7 -ray narrow line Seyfertl are in RPD, while BE Facs 
are basically in the GPD regime. We also find that the FS¬ 
RQs and BE Facs are placed along a line from low power/low 
accretion to high power/high accretion. This is well-known 
“blazar main sequence”, where FSRQs have a strong disc and 
evolve to poorly accreting BE Facs (Cavaliere & D’Elia 2002; 
Bdttcher & Dermer 2002; Maraschi & Tavecchio 2003; Fos¬ 
chini 2011). We also note that some FBFs are in the mass- 
dependent region (left panel), which suggests that FSRQs and 
FBFs may have relationship. Fi et al. (2010) have studied the 
relation between broadband spectral indices Oox and Q!x 7 for 
Fermi blazars. They found that FSRQs and FBFs occupy the 
same region, which suggests that they have similar spectral 
properties. 

3.2. Ratio of jet power to BZ luminosity versus Eddington 

ratio 

We calculate the BZ luminosity by equation (3) and (4) for 
FSRQs/FRII-HEG and BE Facs/FRI-FEG, respectively. Fig¬ 
ure 4 shows the ratio of jet power to BZ luminosity is as a 
function of Eddington ratio. The maximum BZ luminosity of 
FSRQs/FRII-HEG and BE Facs/FRI-FEG can be calculated 
by equation (3) and (4) with j=l, respectively. In Figure 4, we 
find that most of FSRQs/FRII-HEG and BE Facs/FRI-FEG 
have Pjet > Lbz^’ which may indicate that the BZ mecha¬ 
nism is insufficient to explain the jet power of these objects. 
Foschini (2011) suggested that the Blandford-Znajek mecha¬ 
nism is sufficient to explain the power of BE Facs. However, 
we find that our results are contrary to them. We notice that 
they only used 9 BE Facs in their sample, while we use a large 
sample. At the same time, we also find that they used the up¬ 
per limits on the broad line luminosity to the BE Facs. These 
differences may lead to the different results. The above results 
may suggest that there is need to invoke alternative or addi¬ 
tional mechanism to explain the jet power of these objects. 
A hybrid model (i.e., BPh-BZ mechanism) has been recently 
adopted by Garofalo et al. (2010) to explain the observed dif¬ 
ferences in the AGN with relativistic jets. We also find a small 
number of FSRQs and BE Facs have Pjet < Lg“, which may 
indicate that the BZ mechanism is sufficient to explain the jet 
power of this small number of FSRQs and BE Facs. 

3.3. The FRI/FRIIDichotomy 

Meier (1999) given the dividing line between FRI and FRII 
sources of Fedlow and Owen (1996), which can be expressed 
as 

logPrad = -0.66Mr -I-10.35 (5) 

where P^ad is the observed radio power at 1.4 GHz in units 
of W Hz“^ and Mr is the absolute optical R-band magnitude 
of the host galaxy. Mclure and Dunlop (2002) have gotten an 
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Figure 3. Jet power versus black hole mass (left panel) and jet power versus Eddington ratio (right panel), respectively. The 3cr confidence band of the fits is 
shown as the grey shaded regions. The grey shaded regions indicate only the FSRQs and FRII-HEG in the logPjet — logM plane; the grey shaded regions 
indicate all sources in the log Pjet — log Lboi/LEdd plane. The black line is the best least square fits. 



Figure 4. Ratio between the jet power and the calculated luminosity accord¬ 
ing to the Blandford-Znajek theory. The dotted line is j=l. The meanings of 
different symbols are as same as Fig.3. 

empirical relation between Mr of the host galaxy and central 
BH mass, which can be expressed as 

log(MBH/M 0 ) = -0.5(±0.2)Mr - 2.96(±0.48) ( 6 ) 

Willott et al. (1999) pointed out that the jet power can be 
estimated by the following formula 

Qjet 3 X 1038f3/2L®/^_^^,W (7) 

where L 151 ext is the extend radio luminosity at 151 MHz in 
units of 10^® W Hz“^sr“^, and the normalization is uncertain 
and introduce the factor f. Blundell and Rawlings (2000) sug¬ 
gested that the factor f is most likely in the range of 10-20. 


Wu and Cao (2008) obtained the dividing line between FRJ 
and FRII radio galaxies in the Qjet — Mrh plane by using 
equation (5)-(7), Xu et al. (2009) have used different cosmo¬ 
logical parameters that were presented in this Letter to modify 
it, which can be expressed as 

log Qjet (ergs“^) = 1.13 log Mbh(M©) 

-f33.18-f 1.51ogf 

In Figure 5, we plot the relation between the black hole 
mass (M) and jet power (Pjet) for FSRQs and BL Lacs. Fig¬ 
ure 5 shows that the FRI/FRJI diving line was given by Led- 
low and Owen (1996) roughly separates the FSRQs from 
BL Lacs in the M — Pjet plane with f= 10/20, which sup¬ 
ports unihcation schemes of AGN. Blundell and Rawlings 
(2000) suggested that f 10 is likely the consequence of the 
evolution of magnetic field strength as radio sources evolve, 
which would mean that the bulk kinetic (jet) and radiative out¬ 
puts of radio loud AGNs are similar magnitude (Willott et al. 
1999). Willott et al. (1999) found that the ionization lumi¬ 
nosity of radio galaxies is roughly equal to the jet power for 

f=20, which corresponds to Lion/LEdd ~ 2.5 x 10“^ for typ¬ 
ical black hole mass Mbh ~ lO^ ^^^ ^M© in sample of Led- 
low and Owen (1996). It is still unclear why the jet power of 
FRII radio galaxies is always above this dividing line, which 
is beyond the scope of this work. From Fig.5, we also hnd 
that most of LBLs are above the dividing line, while most 
of HBLs are below the dividing line, which is similar to the 
FRI/FRII division. The HBLs have relatively lower jet power 
than LBLs and FSRQs. These results may provide the useful 
clues to investigate the relationship between FSRQs, LBLs 
and HBLs. 

In Figure 6, we show that the distributions of Edding¬ 
ton ratios for BL Lacs and FSRQs, which exhibit a bi- 
modal nature. The FSRQs are clearly separated from the 
BL Lacs at Lboi/LEdd 0.01. Most of BL Lacs have 
Lboi/pEdd < 0.01, while all FSRQs have Lboi/pEdd > 0.01. 
These results that the bimodal behavior of the distribution 
may imply different accretion modes in BL Lacs and FS- 
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Figure 5. Jet power as a function of black hole mass for various classes. 
The two solid line represent the Ledlow-Owen dividing line between 
FRI and FRII radio galaxies given by equation (8) with f=20[top] and 
f= 10[bottom],respectively. The meanings of different symbols are as same 
as Fig.3. 



Figure 6. Distributions of Eddington ratios for Fermi FSRQs (solid line) and 
BL Lac objects (dot line). 

RQs. From Figure 6, we also find that the transition of ac¬ 
cretion states happens at Lboi/LEdd ~ 0.01. Narayan and 
Yi (1995) suggested that the ADAF may be presented when 
Lboi/LEdd < 0.01, which suggested that ADAFs are pre¬ 
sented in BL Lacs and standard thin disks are in FSRQs. A 
similar explanation is proposed to explain the FRl/FRll divi¬ 
sion, in which ADAFs would be presented in FRI radio galax¬ 
ies and the standard thin disks are in FRll radio galaxies (e.g., 
Ghisellini & Celotti 2001; Wu & Cao 2008; Xu et al. 2009). 

3.4. The jet power vi broad line luminosity 



Figure 7. Jet power as a function of broad line luminosity for vari¬ 
ous classes. Shaded colored areas coiTespond to 1, 2 and 3 cr (verti¬ 
cal) dispersion, a = 0.5 dex. The black line is best least square fits 
(logPjet = 0.45 log Lblr + 25). The meanings of different symbols are 
as same as Fig.3. 

The properties of relativistic jets are thought to be closely 
linked with the properties of both the accretion disk and the 
black hole in AGNs. However, the origin and formation of 
relativistic jets are still an unsolved mystery in astrophysics. 
In current theoretical models of the formation of jet, power 
is generated via accretion and extraction of rotational energy 
of disc/black hole (Blandford & Znajek 1977; Blandford & 
Payne 1982), and then converted into the kinetic power of the 
jet. The relation between jet power and broad line luminosity 
is shown in Figure 7. We use a linear regression to analyze 
the relationship between jet power and broad line luminosity. 
From figure 7, we find a significant correlation between them 
(r=0.7081, P<0.0001). This result supports the close connec¬ 
tion between jet power and accretion. The result is consistent 
with other authors (e.g., Rawlings & Saunders 1991; Falcke & 
Biermann 1995; Cao & Jiang 1999; Sbarrato et al. 2012,2014; 
Xiong & Zhang 2014). According to SED modeling, Ghis¬ 
ellini et al. (2014) also found a close connection between jet 
power and accretion for Fermi blazars. 

4. CONCLUSIONS 

In this work, we have explored the relationship between FS¬ 
RQs, BL Lacs, FRI-LEG and ERII-HEG radio galaxies. Our 
main results are the following: 

(i) According to study the Qjet - M and Qjet - Lboi/LEdd 
plane, we find two different regimes. One is the jet power de¬ 
pends on the accretion and another is the jet power depends 
on the black hole mass. ESRQs and ERII-HEG are in the RPD 
regime (mass-dependent regime), while BL Lacs and ERI- 
LEG are in the GPD regime (accretion-dependent regime). 

(ii) Through studying the ratio of between the jet power 
and the calculated the luminosity according to the Blandford- 
Znajek theory versus the Eddington ratios plane, we find that 
the most of both ESRQs and ERII-HEG and both BL Lacs 
and ERI-LEG have Pjet > Lg'l’', which suggests that the BZ 
mechanism is insufficient to explain the jet power of these ob¬ 
jects. There is need to invoke alternative or additional mech- 
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anism to explain the jet power of these objects. We also find 
a small number of FSRQs and BL Lacs have Pjet < 
which may indicate that the BZ mechanism is sufficient to ex¬ 
plain the jet power of this small number of FSRQs and BL 
Lacs. 

(iii) We find that the FSRQs are roughly separated from the 
BL Lacs by the Ledlow-Owen FRI/FRII dividing line in the 
M — Pjet plane. This result supports the unification schemes 
of AGNs. 

(iv) The Eddington ratios Lboi/LEdd of BL Lacs are sys¬ 
tematically lower than those of FSRQs in our sample with 
a clearly division at Lboi/LEdd ~ 0.01. The Eddington dis¬ 
tributions of BL Lacs and ESRQs exhibit a bimodal nature, 
which imply that the accretion mode of BL Lacs may be dif¬ 
ferent from that of ESRQs. 

(v) We find a significant correlation between jet power and 
broad line luminosity in AGNs, which supports the close con¬ 
nection between jet power and accretion. 
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